TECHNICAL  NOTE 

Effect  of  Suspended  Sediment  on 
Acoustic  Detection  Using  Reverberation 


ABSTRACT 


Sonar  operates  by  ensonifying  a  broad  swath  of  the  seabed  using  a  line  array  of 
acoustic  projectors  with  acoustic  backscattering  from  the  ensonified  sediment.  The  sus¬ 
pended  sediment  layer  affects  the  sonar  imagery  through  the  volume  scattering  strength. 
Understanding  the  acoustic  characteristics  of  the  suspended  sediment  layer  can  aid  the 
Navy  in  detecting  sea  mines  with  sonar  imagery.  In  this  study,  the  Navy’s  Comprehen¬ 
sive  Acoustic  Simulation  System  is  used  to  investigate  such  an  effect.  A  range  of  criti¬ 
cal  values  of  volume  scattering  strength  for  buried  object  detection  is  found  through 
repeated  model  simulations. 
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INTRODUCTION 

coustic  detection  of  undersea  objects 
is  difficult  due  to  the  uncertain  en¬ 
vironment  (Chu  et  al.,  2002,  2004)  and 
even  more  difficult  when  the  objects  are 
buried  in  the  seabed.  First,  sediments  gen¬ 
erate  high  backscattering  noise  due  to  het¬ 
erogeneous  scatters  within  the  sediments 
clouding  the  object. 

Second,  the  acoustic  wave  attenuation 
in  sediments  is  much  higher  than  in  water. 
Acoustic  shadows  make  the  buried  targets 
absent  in  the  sonar  images  due  to  diffrac¬ 
tion  around  the  target,  transmission  through 
the  target  and  relatively  high  acoustic  noise 
due  to  backscattering  from  sediments  sur¬ 
rounding  the  target.  Classification  of  bur¬ 
ied  targets  is  also  more  difficult  since  there 
are  no  shadows,  and  the  images  do  not  con¬ 
tain  much  information  about  target  shape 
since  scattering  from  oblique  target  surfaces 
is  not  detectable. 

Acoustic  images  of  buried  targets  pri¬ 
marily  consist  of  echoes  from  the  target  sur¬ 
faces  that  are  normal  to  the  incident  acous¬ 
tic  ray  path.  Target  surfaces  with  an  oblique 
aspect  to  the  incident  ray  path  will  back- 
scatter  much  less  energy  at  the  lower  oper¬ 
ating  frequencies  of  sub-bottom  profilers 
since  the  acoustic  wavelength  is  much 
longer  than  the  surface  roughness  of  most 
targets  of  interest. 


The  suspended  sediment  layer  occupies 
the  lower  water  column.  Presence  of  the  sus¬ 
pended  sediments  creates  a  volume  scatter¬ 
ing  layer  which  affects  acoustic  detection. 
Understanding  the  acoustic  effects  of  the 
suspended  sediment  layer  leads  to  the  de¬ 
velopment  of  acoustic  sensors  with  capabil¬ 
ity  to  scan  the  seafloor  and  to  detect  ord¬ 
nance  such  as  sea  mines. 

2.  Comprehensive  Acoustic 
Simulation  System 

Sonar  equations  provide  guidelines  for 
system  design  (Urick,  1983).  The  govern¬ 
ing  equation  for  beam  patterns  with  domi¬ 
nating  volume  reverberation  is  given  by 

SL-TL-TL  +TS-RL  =  SNR 
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(1) 

where  SL  is  the  source  level;  TL.  is  the 

t 

transmission  loss  of  the  incident  wave; 
TL  is  the  transmission  loss  of  the  re- 

r 

fleeted  echo;  TS  is  target  strength;  RL  is 
the  reverberation  level  of  sediments;  SNR 
is  the  signal-to-noise  ratio  of  the  sonar 
data.  The  transmission  losses  of  the  inci¬ 
dent  and  reflected  waves,  TL.  and  TLr, 
account  for  spherical  loss  (such  as  spheri¬ 
cal  spreading),  acoustic  attenuation,  and 
boundary  loss. 


Volume  reverberation  depends  on  the 
physical  properties  of  the  water  column. 
With  the  presence  of  a  suspended  sediment 
layer,  large  quantities  of  sediment  remain  in 
the  water  column  and  significantly  affect  the 
acoustic  transmission  in  the  water.  The 
denser  the  suspended  sediment  layer  is,  the 
harder  it  would  be  for  sonar  to  penetrate 
through  the  water.  Moreover,  suspended 
bottom  sediment  layer  increases  the  density 
of  the  lower  water  column  and  in  turn 
changes  the  sound  velocity  profile  and  pre¬ 
vents  acoustic  energy  from  reaching  a  pos¬ 
sible  buried  object. 

Reverberation  can  be  used  to  represent  the 
effects  of  a  suspended  sediment  layer  on  acous¬ 
tic  detection.  The  differential  form  of  the  re¬ 
verberation  can  be  obtained  from  integration 
over  the  ensonified  area  (Keenan,  2000), 

d (RL)  =  SL  +  10  log  (&4)  +  55  -  7Z  -  TLf 
+  BPt  +  BPr  ‘  (2) 

where  SA  is  the  scattering  area;  SS  is  the  scat¬ 
tering  strength  per  unit  area;  TL  is  the  trans¬ 
mission  loss  to  scatterer;  7X  is  the  transmis¬ 
sion  loss  from  scatterer;  BP  is  the  beam 

t 

pattern  to  scatterer;  and  BP^ is  the  beam  pat¬ 
tern  from  scatterer.  The  most  important 
design  criterion  for  detecting  mine-like  ob¬ 
jects  is  to  maximize  the  SNR,  the  target  echo 
to  scattering  noise  ratio  in  decibels. 
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The  Comprehensive  Acoustic  Simula¬ 
tion  System  (CASS)  is  the  Navy’s  standard 
model  for  acoustic  and  sonar  analysis.  It  in¬ 
corporates  the  Gaussian  Ray  Bundle 
(GRAB)  eigenray  modes  to  predict  range- 
dependent  acoustic  propagation  in  the  150 
Hz  to  100  kHz  frequency  band  (Keenan  et 
al.,  1996;  Keenan  and  Weinberg,  2001). 

CASS  contains  several  equations  for 
sound  speed  conversions.  The  current  OAML 
approved  Sound  Speed  Algorithm  (Chen  and 
Millero,  1977;  Millero  and  Li,  1994)  has  been 
incorporated  into  CASS.  The  Chen-Millero- 
Li  equations  compute  sound  speed  based  on 
depth,  temperature,  and  salinity.  The  Chen- 
Millero-Li,  Wilson  (1969),  and  Leroy  (1969) 
equations  are  all  very  close  in  the  salinity  range 
30  ppt  to  40  ppt.  For  lower  salinities  the 
Chen-Millero-Li  equation  should  be  used. 
Near  shore  off  the  coast  of  Louisiana  may  have 
salinity  variability  especially  near  the  Missis¬ 
sippi  Delta.  Thus,  the  Chen-Millero-Li  equa¬ 
tions  are  used  in  this  study. 

CASS  simulates  the  sonar  performance 
reasonably  well  in  the  littoral  zone  with  given 
environmental  input  data,  such  as  bottom 
type,  sound  speed  profile  and  wind  speed 
and  accurate  tilt  angle  of  the  sound  source 
(Chu  et  al.,  2002, 2005).  CASS  successfully 
modeled  torpedo  acoustic  performance  in 
shallow  water  exercises  off  the  coast  of  South¬ 
ern  California  and  Cape  Cod.  Recently, 
CASS  was  used  to  simulate  mine  warfare 
systems  performance  in  the  fleet  (Keenan  et 
al.,  1996),  and  for  AN/SQQ-32  mine  hunt¬ 
ing  detection  and  classification  sonar. 

CASS  calculates  the  reverberation  in 
nested  do  loops,  seven  deep.  Reverberation  is 
a  function  of  time  and  the  inner  loop  col¬ 
lected  all  the  reverberation  contributions  over 
the  user-requested  sampling  times.  There  are 
two  loops  on  eigenray  paths,  one  for  the  paths 
connecting  the  transmitter  to  the  scattering 
cell  and  the  other  for  the  paths  connecting 
the  scattering  cell  to  the  receiver.  Since  the 
reverberation  is  calculated  in  the  time  domain 
and  there  may  be  contributions  in  the  same 
time  bin  from  different  ranges,  the  next  loop 
increments  the  range.  CASS  combines  the 
five  possible  eigenray  paths  at  each  range  step 
and  decides  if  the  ray  paths  contribute  to  the 
reverberation  time  bin  (Keenan,  2000). 


Test  rays  are  sorted  into  families  of  com¬ 
parable  numbers  of  turning  points  and 
boundary  interactions.  Ray  properties  are 
then  power  averaged  for  each  ray  family  to 
produce  a  representative  eigenray  of  that 
family.  Target  echo  level  and  reverberation 
level  are  computed  separately,  and  sub¬ 
tracted  to  get  the  signal-noise  ratio  in  the 
absence  of  additive  ambient  noise — noise 
level  is  typically  power  summed  with  the 
reverberation  level  for  total  interference.  A 
detection  threshold  is  applied  to  compute 
SE,  and  then  the  peak  signal  is  used  to  de¬ 
termine  SNR  (Figure  1). 

3.  Environmental  and 
Acoustic  Parameters 

A  nearshore  location  with  silty  clay  bot¬ 
tom  on  the  Louisiana  shelf  is  selected. 
Horizontal  extension  of  the  area  is  50  m  x  60 
m.  Total  depth  including  water  and  sedi¬ 
ment  is  around  100  m.  A  hydrographic 
survey  was  conducted  in  the  area 
(Cornelius,  2004).  The  sound  speed  in¬ 
creases  a  little  from  1520  m  s'1  near  the 
ocean  surface  to  35  m  depth,  and  reduces 
drastically  to  1510.5  m  s'1  at  55  m  depth. 
Below  55  m  depth,  the  sound  speed  de¬ 
creases  slightly  with  depth  and  reaches  1510 
m  s'1  at  100  m. 


Suppose  the  sonar  is  towed  at  a  depth  of 
30.4  m  with  a  source  level  of 240  dB.  The  wa¬ 
ter  depth  in  the  area  varies  from  77  to  95  m. 
The  extent  of  the  image  is  approximately  60  m 
in  the  jy-direction  and  50  m  in  the  x-direction. 
The  grain  size  index  for  a  silty  clay  bottom  is  8 
from  the  Naval  Oceanographic  Offices  stan¬ 
dard.  This  parameter  is  related  to  specific 
geoacoustic  parameters  of  bottom  reflection 
(APL-UW  1994).  The  bottom  reflection  ef¬ 
fects  are  modeled  using  the  Rayleigh  model. 

The  water  column  is  assumed  to  be  rela¬ 
tively  clear  above  77  m  depth  with  a  vol¬ 
ume  scattering  strength  of  -95  dB.  The 
source  of  the  water  column  scattering 
strength  is  extracted  from  a  volume  scatter¬ 
ing  strength  database  (CNMOC,  2004).  A 
suspended  sediment  layer  is  present  below 
77  m  depth  characterized  by  a  different  scat¬ 
tering  strength  (-65  dB).  The  source  level  is 
240  dB.  The  sonar  frequency  is  100  kHz. 
The  pulse  length  is  0.00 1  seconds.  The  time 
increment  for  modeling  should  not  exceed 
one  half  of  the  pulse  length  to  achieve  proper 
resolution  of  each  time  step.  Since  the  total 
distance  traveled  from  the  sonar  to  the  end 
of  the  image  is  approximately  50  m,  the  to¬ 
tal  reverberation  time  is  only  0.12  seconds. 
The  maximum  number  of  bottom  and  sur¬ 
face  reflections  is  set  at  30  to  allow  interfer¬ 
ence  with  reflected  eigenrays. 


FIGURE  1 

Steps  taken  to  create  a  total  reverberation  image  from  CASS,  compared  to  the  side  scan  sonar  image. 
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Test  rays  are  sorted  into  families  of  com¬ 
parable  numbers  of  turning  points  and  bound¬ 
ary  interactions.  Target  echo  level  and  rever¬ 
beration  level  are  typically  computed  separately, 
and  then  subtracted  to  get  the  signal-noise  ra¬ 
tio.  In  the  absence  of  additive  ambient  noise, 
the  signal-noise  ratio  is  typically  power  summed 
with  the  reverberation  level  to  calculate  total 
interference.  A  detection  threshold  is  applied 
to  compute  SE,  and  then  the  peak  signal  is 
used  to  determine  SNR  (Figure  1). 


4.  Mine-Like  Object 

Let  a  hollow  mine-like  steel  object  (8  m 
x  5  m  x  2  m)  be  placed  near  the  center  of 
the  area  with  a  height  of  2  m.  For  the  object, 
the  grain  size  index  is  changed  to  -9  (clay) 
and  the  target  strength  is  set  as  -35  dB.  The 
bathymetry  is  also  changed  to  represent  the 
existence  of  the  object.  The  water  depth  is 
kept  the  same  (87  m)  in  the  vicinity  of  the 
object,  and  changed  into  85  m  over  the  object 
(Figure  2). 

CASS  is  integrated  with  the  sonar  param¬ 
eters,  sound  speed  profile  (SSP),  bottom  type, 
bathymetry,  and  the  scattering  characteristics. 
The  model  output  of  the  seafloor  reverbera¬ 
tion  is  used  to  represent  the  model  generated 
sonar  imagery  (MGSI).  An  increase  of  the 
volume  scattering  strength  in  the  lower  water 
column  reflects  the  presence  of  the  suspended 
sediment  layer.  Clearly,  the  object  is  visible  in 
the  reverberation  imagery.  Since  MGSI  (Fig¬ 
ure  3)  is  a  replica  of  the  sonar  image,  the  ef¬ 
fect  of  suspended  sediment  on  acoustic  de¬ 
tection  may  be  studied  using  MGSI. 


5.  Effect  of  Suspended 
Sediment 

Suspended  sediment  increases  the  vol¬ 
ume  scattering  strength.  To  simulate  its  ef¬ 
fect,  the  volume  scattering  strength  is  in¬ 
creased  by  an  increment  of  5  dB  from  the 
value  of  -65  dB  below  78  m  depth  while 
keeping  the  volume  scattering  strength  con¬ 
stant  (-95  dB)  above  78  m  depth,  and  CASS 
is  integrated  with  increasing  the  value  of  the 
volume  scattering  strength  to  investigate  the 
effect  of  the  suspended  sediment  on  the 
object  detection. 


FIGURE  2 

Bathymetry  with  mine-like  object. 


Depth  (meters)  with  mine-like  object 


FIGURE  3 

Bottom  reverberation  with  mine-like  object  inserted.  Here,  the  horizontal-axis  is  cross  track  indices  represented 
by  time  and  the  vertical-axis  is  along  track  indices  represented  by  distance. 
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The  following  procedure  is  used  to  de¬ 
termine  the  criterion  of  the  volume  scatter¬ 
ing  strength  for  the  object  absent  from 
MGSI.  If  the  simulated  mine-like  object  is 
still  visible  in  MGSI  after  increasing  the  vol¬ 
ume  scattering,  the  simulated  suspended 
sediment  layer  is  not  strong  enough  to  rep¬ 
resent  a  layer  that  would  prevent  a  mine¬ 
like  object  from  sonar  detection.  So  the  vol¬ 
ume  scattering  strength  is  increased  further. 
This  procedure  is  performed  until  the  mine¬ 
like  object  is  no  longer  visible.  During  the 
simulation,  the  water  column  is  assumed  to 
be  relatively  clear  above  77  m  depth,  with 
an  initial  suspended  sediment  layer  charac¬ 
terized  by  a  slightly  stronger  scattering 
strength  below  78  m. 

Sediment  in  the  water  column  would 
increase  the  volume  scattering  and  ulti¬ 
mately  the  volume  reverberation.  For 
MGSI  with  an  object,  a  critical  value  of 
the  volume  scattering  strength  can  be  de¬ 
termined  below  78  m  depth  to  render  the 
mine-like  object  undetectable.  Volume  at¬ 
tenuation  and  changes  in  the  sound  veloc¬ 
ity  profile  will  also  have  an  effect,  but  they 
are  not  addressed  in  this  work. 

As  the  volume  scattering  strength  of  the 
sediment  layer  (below  78  m  depth)  in¬ 
creases  to  a  value  of  -30  dB,  the  object  be¬ 
comes  nearly  undetectable.  The  CASS 
modeling  continues  with  the  increase  of  a 
smaller  increment  of  1  dB.  The  mine-like 
object  is  completely  obscured  (Figure  4)  by 
the  suspended  sediment  layer  at  a  value  of 
-22  dB  (below  78  m  depth),  which  is  taken 
as  the  threshold. 

This  threshold  (-22  dB)  is  large  com¬ 
pared  to  existing  observational  data.  For  ex¬ 
ample,  Kringel  et  al.  (2003)  measure  the 
acoustic  volume  scattering  strength  in  West 
Sound,  Orcas  Island,  Washington,  using  the 
benthic  acoustic  monitoring  system.  It  is  a 
bottom-mounted,  radially  scanning  sonar 
designed  to  record  high-frequency  scatter¬ 
ing  from  the  seafloor.  The  acoustic  volume 
scattering  strength  is  measured  from  22  m 
of  water  column  with  0.5  m  vertical  and  2 
min  temporal  resolution.  Their  measure¬ 
ment  shows  that  the  volume  backscattering 
strength  varies  between  -75  dB  to  -25  dB. 
This  indicates  that  the  threshold  (-22  dB)  is 
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in  the  high  end,  which  implies  that  high 
suspended  sediment  density  is  needed  to 
completely  block  the  mine-like  object  from 
acoustic  detection. 

6.  Conclusions 

(1)  CASS  is  used  to  investigate  the  effect  of 
suspended  sediment  on  detecting  a  mine¬ 
like  object  in  the  silty  clay  bottom  at  the 
Louisiana  shelf  with  water  depth  around  1 00 
m.  Hydrographic  and  meteorological  sur¬ 
veys  were  conducted.  The  environmental 
data  (wind  and  SSP)  are  taken  as  the  input 
into  CASS  to  simulate  the  sonar  image  with 
a  mine-like  object  present  through  its  rever¬ 
beration  characteristics. 

(2)  A  threshold  of  volume  scattering  strength 
(-22  dB)  for  the  mine-like  object  detection 
is  found  through  repeated  model  simula¬ 
tions.  When  the  volume  scattering  strength 
increases  to  the  threshold,  the  mine-like 
object  is  acoustically  undetectable.  It  is  noted 


that  this  threshold  is  very  large  and  valid  only 
for  this  case  and  that  the  purpose  of  this 
study  is  to  show  the  methodology  rather  than 
to  provide  an  accurate  threshold. 

(3)  While  CASS  is  useful,  several  short¬ 
falls  remain  in  this  study.  First,  the  pro¬ 
cess  by  which  the  environment  and  the 
object  are  modeled  is  cumbersome.  Sec¬ 
ond,  the  appropriate  volume  scattering 
strength  for  the  buried  object  should  be 
given.  A  thorough  study  is  suggested  on 
the  relationship  between  the  suspended 
sediment  layer  density  and  type  (e.g.,  sand, 
silt  or  clay),  particle  density  in  the  layer, 
associated  volume  scattering  strength  and 
attenuation,  and  changes  in  the  sound 
speed  profile. 
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FIGURE  4 

Reverberation  plot  of  bottom  with  mine  object  inserted.  The  horizontal  axis  is  cross  track  indices  represented  by 
time  and  the  vertical  axis  is  along  track  indices  represented  by  distance  with  mine  object  inserted  and  sus¬ 
pended  sediment  layer  with  -22  dB  volume  scattering  strength. 
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